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ABSTRACT 

Scaling relations of observed galaxy cluster properties are useful tools for constraining cosmological 
parameters as well as cluster formation histories. One of the key cosmological parameters, as, is 
constrained using observed clusters of galaxies, although current estimates of as from the scaling 
relations of dynamically relaxed galaxy clusters are limited by the large scatter in the observed cluster 
mass-temperature (M — T) relation. With a sample of eight strong lensing clusters at 0.3 < z < 0.8, 
we find that the observed cluster concentration-mass relation can be used to reduce the M — T scatter 
by a factor of 6. Typically only relaxed clusters are used to estimate as, but combining the cluster 
concentration- mass relation with the M — T relation enables the inclusion of unrelaxed clusters as well. 
Thus, the resultant gains in the accuracy of as measurements from clusters are twofold: the errors on 
as are reduced and the cluster sample size is increased. Therefore, the statistics on as determination 
from clusters are greatly improved by the inclusion of unrelaxed clusters. Exploring cluster scaling 
relations further, we find that the correlation between brightest cluster galaxy (BCG) luminosity and 
cluster mass offers insight into the assembly histories of clusters. We find preliminary evidence for 
a steeper BCG luminosity - cluster mass relation for strong lensing clusters than the general cluster 
population, hinting that strong lensing clusters may have had more active merging histories. 
Subject headings: cosmological parameters - clusters: individual (3C 220, A 370, CI 0024, CI 0939, 
CI 2244, MS 0451, MS 1137, MS 2137) - dark matter - galaxies: evolution - 
galaxies: formation - gravitational lensing 
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1. INTRODUCTION 

As the most massive bound systems known, galaxy 
clusters provide an important link in understanding the 
composition and growth of structure in the Universe. 
Clusters follow a variety of observational scalings of mass 
with temperature, luminosity, or cluster counts, and 
these scalings are sensitive to cosmological parameters 
including the matter density parameter f2 m , the cosmo- 
logical constant density parameter the dark energy 
equation-of-state parameter w, and the normalization 
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1091) . Such constraints from galaxy clusters comple- 
ment the constraints on cosmological parameters from 
Type la supernovae and cosmic microwave background 
observations. 

However, useful galaxy cluster constraints on cosmo- 
logical parameters depend primarily on accurate deter- 
minations of cluster masses. Observationally, cluster 
masses are typically measured in one of three ways. 

A long-established method for determining cluster 
masses employs the virial theorem and the measure- 
ment of velocities of the galaxies that constitute the clus- 
ter. Based on the three assumptions that the cluster is 
in virial equilibrium, the galaxy distribution efficiently 
traces the cluster mass distribution, and the velocity dis- 
persions a of the galaxies are isotropic, the cluster mass 
contained within a radius r is estimated M ~ a 2 r/G. 
However, these mass estimates may be biased as a re- 



sult of galaxy velocity anisotropics or if the galaxy dis- 
tribution do es not follow the total mass distribution (e.g., 
lBailevlll982h . 

A second method uses cluster X-ray emission as a 
tracer of cluster masses. The hot intracluster gas, which 
is the dominant baryonic component of a cluster and is 
typically twice the mass of the total mass of the galaxies 
in a cluster, emits X-rays via bremsstrahlung radiation 
and atomic line emission. With the temperature T and 
radial density p(r) profiles determined from X-ray spec- 
tra and surface brightness distributions, the cluster mass 
is given by M ~ r 2 /p(r) d{—pT)/dr. This method as- 
sumes that the intracluster gas is s pherically distribute d 
and is in hydrostatic equilibrium (|Evrard et al.l Il996h . 
However, these assumptions may be incorrect. If the gas 
distribution is not spherical, X-ray mass estimates will 
be biased by projection effects. Many galaxy clusters are 
also not in hydrostatic equilibrium, in particular dynam- 
ically unrelaxed clusters that are undergoing mergers. 
There is evidence that the bias of hydrostatic equilib- 
rium mass is linked to the dynamical state of the galax y 
cluster (e.g., lAndersson et aJJ 12009b IZhang et all I2009Q . 
In addition, the hot gas of galaxy clusters with buoyant 
bubbles near their cores migh t indicate a departure from 
hydrostatic equilibrium (e.g.. IChurazov et al.ll2001l ). 

The most direct estimates of cluster masses employ 
gravitational lensing distortions of background galaxies. 
This technique is free of assumptions about the dynam- 
ical state of the cluster, which enables it in principle to 
yield more consistent mass estimates, though it is also 
sensitive to projection effects. More accurate cluster 
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mass estimates can in turn provide tighter constraints 
on cosmological parameters, and therefore it is of key im- 
portance to reduce the errors in cluster mass estimates. 

For example, the primary source of error in cluster- 
based determinations of as is the error in the 
mass-temperature re l ation for rel a xed cluster s (e.g., 
iPierpaoli et all 12001 iHenrvl [200l IVoit] f2005h . Re- 
cent studies show that an X-ray independent mass ap- 
proach such as gravitational lensing provides a unique 
tool to calibrate t h e mass - tempe r ature r elation (e.g.. 
Smith et all 120051: iMahdavi et all 120081: IZhang et alj 
2008). Here, we use strong gravitational lensing mass 
measurements of a sample of eight strong lensing clus- 
ters at 0.3 < z < 0.8 to accurately measure the galaxy 
cluster mass-temperature relation. We also include the 
effects of cluster concentrations in an effort to further 
reduce the scatter in the cluster mass-temperature rela- 
tion, which would ultimately enable tighter constraints 
on erg. 

In addition to the correlations that exist between clus- 
ter properties, some observational properties of brightest 
cluster galaxies (BCGs) also scale with properties of the 
host clusters. Whereas scalings between cluster proper- 
ties are sensitive to cosmological parameters, scalings be- 
tween BCGs and their host clusters provide constraints 
on BCG formation and the evolution of clusters. 

BCGs are a unique population: they are the most mas- 
sive and luminous galaxies in the Universe. They are 
typically located near the centers of clusters, which sug- 
gests that a BCG's formation history is intricately linked 
to the formation of the cluster itself. However, the for- 
mation of BCGs is still poorly understood. 

BCGs may form after their host clusters assemble in 
one of two ways. First, a BCG may be the first galaxy to 
be dragged in by dynamical friction to the center of the 
dark matter halo destined to become a cluster, where 
it then grows through galactic cannibalism by merging 
with subsequent galaxies t hat fall to the center (e.g . , 
lOstriker fc Tremainel Il975t lHausman fc Ostrikerl [l97l . 
However, this scenario typically requires more than a 
Hubble time to form a BCG because much of the mass 
of the infalling galaxy is tidally stripped, which reduces 
the d ynamical friction effect and slows the infall (jMerrittl 
[19851) . 

BCG formation may also occur after cluster forma- 
tion if the host cluster's central cooling flow forms 
stars at the cluster center a nd those stars build the 
BCG (jCowie fc Binnevl fl977t ). There are several in- 
stances of ongoing or recent star format ion in BCGs 
that occupy cooling-flow c l usters (e.g.. ICardiel et al.l 
19981: iCrawford et al.1 Il999t iHicks fc Mushotzkvl 120051 : 
McNamara et a l. 2006), but it is unclear whether the star 
formation is fueled by the cooling flows o r by cold gas 
broug ht in through recent galaxy mergers ( Bil dfell et all 

In another scenario, BCGs might form in concert with 
their host clusters. A BCG may begin with several galax- 
ies merging together in a group to form a large galaxy, 
and then when groups merge as hierarchical structure 
formation continues, this large galaxy eventually be- 
comes a BCG in a massive cluster fe ^Merrittlrl985l : 
lDubinskil[l998T: iBovlan-Kolchin et al.ll2006t ). 

Here, we examine the correlation between BCG lumi- 
nosity and cluster mass in eight strong lensing clusters at 



0.3 < z < 0.8. This will enable constraints not only on 
BCG and cluster formation in general, but also on how 
the BCGs in strong lensing clusters may have formed 
and evolved differently than BCGs in the general cluster 
population. 

The rest of this paper is organized as follows. In Sec- 
tion [5] we describe the selection of our cluster sample, 
and Section [3] gives the masses, dynamical states, and X- 
ray temperatures for these clusters. In Section |4] we find 
the M — T relation for the relaxed clusters in our sam- 
ple and show how the inclusion of cluster concentrations 
both significantly reduces the scatter in the M — T rela- 
tion and lifts the restriction on cluster dynamical state. 
In Section [5] we identify the BCGs in our sample and 
measure their luminosities. We use these luminosities 
in Section [6] to measure the correlation between BCG 
luminosity and cluster mass, and we find preliminary ev- 
idence that strong lensing clusters may have more ac- 
tive merging histories than the general cluster popula- 
tion. Section [7] presents our conclusions. Throughout 
this paper, we adopt a spatially flat cosmological model 
dominated by cold dark matter and a cosmological con- 
stant (On = 0.3, 0\ = 0.7, h = 0.7). 

2. SAMPLE SELECTION 

We base our sample o n 10 well-known strong lens- 
ing clusters analyzed in IComerford et al.l ((2006). All 
10 clusters have Hubble Space Telescope (HST) imag- 
ing, which make possible the mass determinations and 
photometry measurements central to this paper. How- 
ever, there are no published arc redshifts for two of the 
clusters, CI 0016+1609 and CI 0054-27, which limits 
the strong lensing determination of their cluster masses 
to the unknown factor D s /D\ Sl the ratio of the angular 
diameter distances to the source and between the lens 
and source. Consequently we remove these two clusters, 
and our sample consists of the remaining eight clusters 
at 0.3 < z < 0.8: C1G 2244-02, Abell 370, 3C 220.1, 
MS 2137.3-2353, MS 0451.6-0305, MS 1137.5+6625, 
CI 0939+4713, and ZwCl 0024+1652. 

3. CLUSTER PROPERTIES 

Strong correlations are found between cluster observ- 
ables, and the resultant scaling relations clearly encapsu- 
late key information about cosmological parameters and 
the assembly history of clusters. Cluster masses are a 
component of many cluster scaling relations, and we mea- 
sure strong lensing masses for our sample of clusters and 
compare these to mass estimates from the distributions 
of cluster X-ray gas. Based on these comparisons and 
other observable properties of the cluster, we determine 
the dynamical state of each cluster as relaxed or unre- 
laxed. We also present cluster X-ray temperatures, which 
are another component of cluster scaling relations. 

3.1. Cluster Strong Lensing Mass Determination 

We model each cluster mass distr ibution with an ellip- 
tical Navarro-Frenk- White (NFW; iNavarro et~ai] I1996L 
Il997h dark matter halo centered on the BCG, using 
the be st-fit NFW parameters found by IComerford et"aLl 
( 2006). Strong lensing arcs with measured redshifts ob- 
served in a cluster con strain its mass distribution, and 
IComerford et~aTl {2006) use the arcs to characterize best- 
fit NFW ellipsoids to each cluster. With the NFW dark 
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matter halos completely defined in this way, we can de- 
termine any cluster radius ta as the radius at which the 
density of the halo is A times the critical density at the 
cluster redshift. 

Lack of information about the clusters' three- 
dimensional shapes prevents us from calculating their 
elliptical masses, but instead we determine the equiv- 
alent mass of a sph erical NFW halo. With the 
iComerford et al.l (|2006l ) best-fit scale convergence k s and 
scale radius r s , we estimate the cluster mass within ra- 
dius ta as 



M A = 47r£ cr itK s r- s 



ln(l + x) 



1 



(1) 



where x = r&/r s and £ C rit is the critical surface mass 
density, defined as 



A 
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which depends on the angular diameter distances A,s,is 
from the observer to the lens, to the source, and from 
the lens to the source, respectively. 

We estimate the errors in mass by propagating the 
errors in the bes t -fit N FW parameters. As detailed in 
IComerford et al.l (|2006T ) these errors are quite small but 
are realistic, because the reproduced lensed image is sen- 
sitive to slight variations in a parameter's value. How- 
ever, we note that these errors are relevant only to the 
choice of lens model and data and do not represent a 
global systematic uncertainty. 

We use the method described here to measure the lens- 
ing cluster masses in Table [3 as well as the cluster masses 
M 2 oo and M 250 o in Table [2] 

3.2. Dynamical State of Clusters: 
Relaxed vs. Unrelaxed 

Since one of our aims is to measure the mass- 
temperature relation for relaxed lensing clusters, we must 
determine which of the eight clusters in our sample are 
dynamically relaxed. X-ray cluster mass estimates are 
based on the assumption that the cluster is in hydro- 
static equilibrium, and if a cluster is relaxed it is also 
in hydrostatic equilibrium. Therefore, X-ray mass mea- 
surements for relaxed clusters should be accurate and 
consistent with lensing mass measurements. 

We use X-ray mass estimates from the literature, where 
the X-ray masses are measured for each cluster at two or 
three different radii. For each cluster, Table [T] gives the 
lensing mass and X-ray mass measured within the two 
or three different cluster radii. Table Q] also shows the 
lensing mass to X-ray mass ratio and the reduced % 2 
of the comparison of lensing and X-ray masses. For six 
clusters, at all radii at which masses were measured, the 
ratio of lensing mass to X-ray mass is consistent with 
unity and the reduced \ 2 1S S 1) suggesting that these 
six clusters could be relaxed. Additional observational 
evidence in § 13.2.11 and § 13.2.21 shows that four of these 
six clusters are relaxed, while the remaining two clusters 
are unrelaxed. 

For at least one of the radii considered, the two clusters 
MS 2137-23 and CI 0939+4713 each exhibit lensing to 
X-ray mass ratios that are inconsistent with unity and 
reduced x 2 that are greater than unity, which is evidence 



that the clusters are unrelaxed. We measure masses for 
MS 2137—23 within three different radii, and within one 
of these radii the mass ratio is inconsistent with unity 
and the reduced \ 2 is greater than unity. However there 
is opposing evidence, given in § 13.2. 1\ that characterizes 
MS 2137-23 as a relaxed cluster. For CI 0939+4713, 
the mass ratios measured at both radii considered are 
inconsistent with unity and both reduced \ 2 are much 
greater than unity, suggesting CI 0939+4713 may be an 
unrelaxed cluster. In § 13.2. 21 we present more evidence in 
support of this conclusion. 

Additional information about the dynamical state of 
a cluster can be found in its X-ray emission map. For 
example, the position of the BCG relative to the peak in 
the cluster's X-ray profile may be evidence of a cluster's 
dynamical state: if the two are coincident the cluster 
is likely relaxed, otherwise it is likely unrelaxed. The 
centroid shift is one means of quantifying this positiona l 
difference (e.g., iMohr et al.l 119931: Ueltema et alll2008[ ). 
Additionally, a smooth distribution of X-ray gas indi- 
cates the cluster is likely in a relaxed state. However, if 
the X-ray gas is distributed irregularly or shows evidence 
of shocks or substructure, the cluster is likely unrelaxed 
and undergoing a merger. Below we examine evidence 
for the dynamical state of each cluster individually and 
label each cluster as relaxed or unrelaxed (these labels 
are also given in Table [1]). We first discuss the four re- 
laxed clusters, then the four unrelaxed clusters. 

3.2.1. Relaxed Clusters 

• CI 2244-02: We find that X-ray and lensing 
mas ses for CI 2244 —02 are consistent (Table [lj 
and lOta et al.l (|1998l ) also find consistent X-ray and 
lensing masses, suggesting that hydrostatic equilib- 
rium is a valid assumption for CI 2244—02 and that 
it is a relaxed cluster. 

• 3C 220.1: The radial profile of X-ray emission 
from 3C 220.1 shows no sign of irregularity and 
the profile is well-fit by a model assuming hydro- 
static equilibrium , which suggest tha t 3C 220.1 is 
a relaxed cluster (fWorrall et al.ll200l1) . 

• MS 2137-23: The X-ray and strong lensing 
masses of M S 2137—23 are in good agreement 
(Allen 1998), indicating that it is in a relaxed state. 
Many relaxed clusters also have cooling flows , such 
as the massive cooling flow in MS 2137-23 (|Allenl 
[19981: IWull2000h . 

• MS 1137+66: The cluster MS 1137+66 not 

only has consistent X-ray and weak lensing masses 
(Table [U), but also has a small centroid shift 



Maughan et a l. 2008) and may host a moder- 
ate cooling flow ([Donahue et al.1 1X9991 ) . In addi- 
tion, Sunyaev Zel'dovich obser vations of the cluste r 
show no obvious substructure (jCotter et al.l l2002'). 
These properites connote that MS 1137+66 is a re- 
laxed cluster. 



3.2.2. Unrelaxed Clusters 

• Abell 370: Abell 370 hosts two cD galaxies, 
and there are X-ray peaks centered on each cD 
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TABLE 1 

Comparisons between strong lensing and X-ray cluster mass estimates. 
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TABLE 2 

Cluster lensing masses and X-ray temperatures. 
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(|Mellier et al.lll994l ). The two cD galaxies are mov- 
ing relative to each other at 1000 km s" 1 , signaling 
that Abell 370 is an unrelax ed cluster undergoing 
a merger dKneib et al.lll993h . 

• CI 0939±4713: X-ray observations of 
CI 0939+4713 show evidence for substructure 
(jSchindler fc Wam bsganss 1996), and the disagree- 
ment between lensing and X-ray masses shown 
in Table [Q further suggests that CI 0939±4713 is 
not in hydrostatic equilibrium. These observations 
indicate CI 0939±4713 is an unrelaxed cluster. 

• CI 0024+17: The two dark matter clumps 
near the center of CI 0024+17 are separated 
in r edshift, implying that it is a merging clus- 
ter (|Nataraian et al.l [20091 . There is additional 
evidence for substructure in CI 0024+17 in its 
mass models, which require substructure to pro- 
duce a good fit to th e cluster's lensing arcs 
(|Broadhurst et al.l |2000D . The redshifts of the 



member galaxies are distributed bimodally, fortify- 
ing the evidence that CI 0024+17 ma y have under- 
gone a merger with another cluster (jCzoske et al.1 
2002). The evidence implies that CI 0024+17 is an 
unrelaxed cluster. 

• MS 0451.6-0305: The distribution of mass 
within the central 1' of MS 0451.6-0305 is not 
smooth, and the centroid shift i ndicates the BCG 
is not located at the X -ray peak (jBorvs et al.ll2004l : 
iMaughan et al.l 120081 ). These observations suggest 
that MS 0451.6-0305 is unrelaxed. 

3.3. Cluster X-ray Temperatures 

The temperature of the intracluster medium is com- 
monly measured using its X-ray emission in one of 
several ways: through fits to the cluster's observed 
X-ray spectrum (yielding the spectroscopic tempera- 
ture T s ), through weighting by the mass of the gas 
element (yielding the mass- weighted temperature T m ), 
or through weighting by the emissivity of the gas 
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element (yielding the emission-weighted temperature 
Tew.)- However, the spectroscopic temperature T s is 
systematically lower than the mass-weighted temper- 
ature T m and the emission-weighted temperatu re T em 
(jMathiesen fc EvrardlfeOOlUMazzotta et al.ll2004| ). so an 
accurate temperature comparison across different clus- 
ters requires consistent temperature measurements. 

To ensure that the cluster temperatures we use for our 
sample are as consistent as possible, we use the mean 
cluster temperatures derived from the single -temperature 
model fits of ASCA data in iHornerl (|2001t ). This large, 
homogeneous catalog of spectroscopic cluster tempera- 
tures includes six of our clusters, and for the remain- 
ing two clusters, 3C 220.1 and CI 0939+4713, we remain 
as consistent as possible by using spectroscopic temper- 
atures from single-temperature fits. Table [5] gives the 
cluster temperatures and the corresponding references. 
We note that none of the temperatures we use apply cor- 
rections for cool cores at the cluster centers. 

Some clusters in our sample also have temperature 
measurements from Chandra and XMM-Newton data. 
Specifically, 3C 220.1 has a Ch andra tempera t ure o f 
8.5±|;| keV within 10 - 45" (IWorrall et al.l l200l : 

MS 0451.6-0305 has a Chandra temp erature of 6.7±g;| 
keV within r 500 (|Maughan et al.ll2008l ): MS 1137.5+6625 
has a Chandra temperature of 5.8^g'g keV within r^Q 
(jMauehan et al.ll2008r i; and CI 0024+1 7 has an average 
Chandra temperature of 4.47j£^ keV (|Ota et al.ll2"00l 
and an X MM- Newton temp erature of 3.52 ± 0.17 keV 
within 3' (| Zhang et al.ll2005f) . Use of Chandra or XMM- 
Newton temperatures could change the results of the 
mass-temperature relation. However, because Chandra 
and XMM-Newton temperatures have been measured for 
only a subset of our sample, and because these temper- 
atures are measured within inconsistent cluster radii, we 
do not use Chandra and XMM-Newton measurements in 
our determination of the cluster mass-temperature rela- 
tion below. 

4. THE MASS-TEMPERATURE RELATION 

Theoretical arguments suggest a correlation between 
cluster mass and X-ray temperature for relaxed clus- 
ters, which provides the link between the gas in a clus- 
ter and its mass. Here we determine the cluster mass- 
temperature relation for relaxed strong lensing clusters, 
and we also explore the correlation between the scatter 
in cluster temperature and the scatter in cluster concen- 
tration to establish a general mass-temperature relation 
that is independent of the dynamical state of the clusters. 

4.1. The M — T Relation for Relaxed Strong Lensing 

Clusters 

A correlation between cluster mass and cluster X-ray 
gas temperature in relaxed clusters is expected as a direct 
consequence of theoretical arguments. If a cluster's X- 
ray gas is in virial and hydrostatic equilibrium, then the 
theoretical expectation is that cluster mass scales with X- 
ray temperature as E(z)Ma = A(A)T 15 , where E(z) = 
H(z)/H = ^fl m (l + z) 3 + n A for a flat Universe, M A 
is the cluster mass within the radius where the mean 
mass density is A times the critical density, and A(A) is 
the A-depcndent normalization. 




CI0939 



10 

kT (keV) 

Fig. 1. — The mass-temperature relation for observed strong 
lensing clusters. Unrelaxed clusters (open circles) are not included 
in the fit, and the relaxed clusters (black points) are fit by a power 
law with slope a = 1.43 (black solid line). The lcr scatter for 
all eight clusters is large, A(log[i?(z)M25oo]) = 0.2 (black dashed 
lines) . Also shown are the other M — T relations for observational 
samples that use spectroscopic temperatures as we do : 17 w eak 
lensing clusters with 3.6 < T s (keV) < 9.8 (Hockstra 20071: red 
dott ed line), 13 relaxed X-ray clusters with 0.7 < T s (keV) < 
8.9 IVikhlinin et al.ll200l blue dash-d otted line), 10 relaxe d X-ray 
clusters with 2.2 < T a (keV) < 8.3 l lArnaud et alj|2005l ; orange 
das hed line), and six rela xed X-ray clusters with 3.7 < T B (keV) < 
8.3 ijArnaud et al.l [20051: green long dashed line). We find that 
our slope is in agreement with both the theoretical expectation 
of a = 1.5 and measurements of a by other observations. For a 
detailed comparison to these and other estimates of the M — T 
relation, see Table [3] 

The critical overdensity A = 2500 is commonly used in 
cluster analyses because in the central regions enclosed 
by ^25oo, Chandra cluster temperature profiles can be 
measured even a t high redshifts (e.g., up to z = 0.9 in 
lAllen et al]l2004t ). The overdensity A = 2500 is therefore 
appropriate for our cluster sample, which extends to z — 
0.8. Using the overdensity A = 2500, we can write the 
cluster mass-temperature relation in power law form as 

Using our sample of four dynamically relaxed lensing 
clusters given in § 13.2.11 a best fit to the power law M — T 
relation yields A = 1.60 ± 3.42 and a = 1.43 ± 1.28, 
consistent with the theoretical expectation of a = 1.5. 
Figure [T] shows this best-fit relation, for which the RMS 
scatter is 360% for all eight clusters and 500% for the 
four unrelaxed clusters. 

We compare with other observations and simulations of 
the M — T relation in Table El including those that used 
spectroscopic temperatures T s and those that used mass- 
weighted temperatures T m . For cases where the temper- 
ature normalization is not 5 keV and/or the mass scaling 
is not 10 14 h^Q Mq , we recalculate A using the published 



6 



Comerford, Moustakas, & Natarajan 



TABLE 3 

Power Law Fits to the M — T Relation. 



A a Method a kT (keV) b Sample Reference 



1.60 ± 3.42 1.43 ± 1.28 SL 4.6 < T B < 6.7 4 relaxed SL clusters 1 

2.0 ±0.29 1. 34±o28 WL 3.6 < T s < 9.8 17 WL clusters 2 

1.79 ±0.07 1.64 ± 0.06 X-ray 0.7 < T B < 8.9 13 relaxed clusters 3 

2.06 ± 0.10 1.58 ± 0.07 X-ray 0.6 < T m < 9.3 13 relaxed clusters 3 

1.69 ±0.05 1.70 ±0.07 X-ray 2.2 < T s < 8.3 10 relaxed clusters 4 

1.79 ± 0.06 1.51 ± 0.11 X-ray 3.7 < T s < 8.3 6 relaxed clusters 4 

1.88 ±0.26 1.52 ±0.36 X-ray 5.6 < T m < 15.3 5 relaxed WL or SL clusters 5 



1.97 ±0.07 1.54 ±0.02 Simulation T m M2500 > 4 X 10 14 hy L M clusters 6 

in hydrodynamics simulation 



References . — ( 1 ) Thi s paper; (2 ) IHoekstral i2007h : (31 1 Vikhlinin et al.l l l200r3l ; f4) IArnaud et all (20051: 
(5) lAllen et all (I2001T) : (6) IKav et all [(20051) . 

a Method used to determine the cluster mass, where SL is strong lensing and WL is weak lensing. 
b Temperature range of the cluster sample, where T s is the spectroscopic temperature and T m is the mass- 
weighted temperature. 




Fig. 2. — The correlation between the difference AT between 
the observed X-ray temperatures and the predicted temperatures 
from the M — T relation and the difference Ac between the mea- 
sured concentrations and the predicted concentrations from the 
c — M relation. The eight strong lensing clusters in our sample are 
represented, and the solid line shows the best-fit line to the data 
AT = (-2.75 keV)Ac - (1.56 keV). The dashed lines show the la 
scatter A (AT) = 0.9 keV. 

slope a, a temperature normalization of 5 keV, and a 
mass normalization of 10 14 h^o M@. To be conservative, 
we assume the fractional error in A is unchanged. 

The observations we compare span varying temper- 
ature ranges, and there is some evidence that the 
M — T relation stee pens for cooler clusters (e.g., 
iNevalainen et al.l l2000HFinoguenov et al.ll2001l ): for ex- 
ample, lArnaud et al.l (j2005f ) find a slope of a = 1.51 for 
clusters with 3.7 < T s < 8.3 keV, which increases to 
a = 1.70 for clusters with 2.2 keV < T s < 8.3 keV. The 
temperature range we probe (4.6 keV < T s < 6.7 keV) is 
likely too small to exhibit a significant change in slope, 
but we lack a large enough statistical sample to test this 



Fig. 3. — The mass-temperature relation, after correcting for 
the scatter in temperature, for observed strong lensing clusters. 
As in Figure [T] open circles represent unrelaxed clusters and black 
points represent relaxed clusters. We adjust the temperature of 
each cluster according to its concentration and the AT — Ac rela- 
tion. The best-fit M — T relation for relaxed clusters, derived in 
§ 14.11 is shown as the solid line. The lcr scatter for all eight clus- 
ters is A(log[i?(2)M25oo]) = 0-1 (black dashed lines), significantly 
smaller than the scatter in the uncorrected M — T relation (see 
Figure [TJ. 



properly. 

We find that our best-fit slope a is consistent with 
both the theoretical expectation and the slopes derived 
by other observations and simulations of clusters. Our 
best-fit normalization A is somewhat lower than, but 
still consistent with, the normalizations found by the 
other observations and simulations. We find that relaxed 
strong lensing clusters follow the same M — T relation as 
relaxed clusters in general. 
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4.2. Correlation between the Temperature Scatter and 
Concentration Scatter 

We have derived anM-T power-law relation for re- 
laxed lensing clusters, but a more general M — T relation 
including both relaxed and unrelaxed clusters may be 
possible if we account for the differences in cluster con- 
centrations. First, we define the virial radius of a cluster 
as the radius r v i r at which the average cluster density 
equals A v ; r (z) times the mean density at the cluster red- 
shift z, where A vir (z) ~ (18tt 2 + 82s - 3 9x 2 )/(l + x) and 
x = £l m (z) — 1 (jBrvan fe Normar]ll998[ ). Using the scale 
radius r s of the best-fit NFW profile to each cluster, the 
cluster concentration is defined as c vlr = r v - lr /r s . 

Since more concentrated clust ers are expec t ed to 
form at higher red shifts (e.g., iNavarro et al.l Il997t 
IWechsler et al.|[2002T) , if the cluster X-ray gas cools with 
time there might be a correlation between high cluster 
concentrations and low cluster temperatures. In addi- 
tion, mergers with other clusters or groups may deplete 
the central mass densities in clusters while shock-heating 
the cluster gas, producing high cluster temperatures for 
low cluster concentrations. Here, we analyze whether 
there is any such correlation between the scatter in tem- 
perature and the scatter in concentration for our sample 
of eight strong lensing clusters. 

Cluster concentrations c vn - and cluster virial masses 
M v j r = M(< r vir ) are determined by strong lensing 
m easurements for each of the cl usters in our sample 
in IComerford fc Natarajanl ()2007l ). The concentration 
c v j r = 16 determined by strong lensing measurements of 
MS 2137.3—2353 is known to be overestimated because 
the cluster's dark matter halo is likely elongated along 
or near the line of sight ()Gavaz zi 2005), so we instead 
use the concentration c V i r = 8.75 der ived from the X- 
ray m ass profile for MS 2137.3-2353 (|Schmidt fe Alien! 
l2007f l. We note that if the lensing concentration were 
used for MS 2137.3-2353, Equation would be AT = 
(-0.07 keV) Ac - (0.49 keV). 

From a sample of 62 galaxy clusters, 

IComerford fc Nataraiar] (|2007l ) find a power-law 
relation between cluster concentration c V i r and cluster 
virial mass M v - lr of 

_ 14.5 ±6.4 / M vir \ -°- 15±0 - 13 

Cvir (1 + z) 1,1.3 x I0^h- X M Q ) ' { > 

where z is the cluster redshift. For each of the eight clus- 
ters in our sample, we calculate the difference Ac between 
the measured concentration and the concentration pre- 
dicted by the above c — M relation. We also calculate the 
difference AT between the measured X-ray temperature 
and the temperature predicted by the M — T relation we 
determined in § 14.11 for the four relaxed clusters. 

Figure [2] shows the results of these AT and Ac calcu- 
lations. The best-fit line to the data is 

AT = (-2.75 keV±0.07 keV)Ac-(1.56 keV±0.49 keV) , 

(5) 

suggesting that indeed higher (lower) temperature clus- 
ters tend to have lower (higher) concentrations. 

4.3. The M — T Relation for All Strong Lensing 
Clusters 

Using the relation between the scatter in cluster tem- 
perature and the scatter in cluster concentration for the 



eight strong lensing clusters (§ I4.2|) . we adjust for the 
apparent dependence of cluster temperatures on clus- 
ter concentrations. We use Ac for each cluster to cal- 
culate its corresponding AT from the best-fit relation 
given in Equation [5j We then subtract this AT from 
the measured temperature to obtain a corrected temper- 
ature T corr , and we illustrate the resultant temperature- 
corrected M — T relation in Figure [3l The figure also 
shows the relation we derived in § 14.11 for the four re- 
laxed clusters, where A = 1.60 and a = 1.43. 

We find that cluster concentration, mass, and X-ray 
temperature are tightly correlated, and as a result in- 
corporating the AT — Ac relation significantly reduces 
the scatter in the M — T relation. Comparing Figure [3] 
to Figure [1] underscores the impact of our temperature 
correction in reducing the scatter in the M — T relation. 
The temperature correction reduces the RMS scatter for 
all eight clusters by a factor of 6, from 360% to 60%, and 
more significantly, reduces the RMS scatter for the four 
unrelaxed clusters by a factor of 30, from 500% to 15%. 
(The RMS scatter for the four relaxed clusters increases 
from 26% to 83%, possibly because the temperatures we 
use do not correct for cool cores at the cluster centers.) 
With the temperature correction, even unrelaxed clusters 
follow the M — T relation we originally derived using only 
the relaxed clusters (§ 14. ip . Therefore, we suggest this 
temperature correction as a tool for establishing a uni- 
versal M — T relation that applies to all galaxy clusters 
regardless of their dynamical state. 

The error in the measurement of as from cluster counts 
depends directly on the error in the cluster M — T re- 
lation; for example, a 25% la uncertainty in the zero 
point of the M — T relation corresp onds to a 10% la 
uncertainty in as (|Evrard et al.l l2002) . Consequently, we 
find that the temperature correction not only reduces 
the scatter in the M — T relation, but also significantly 
reduces the error in the corresponding measurement of 
a s . 

An alternate cluster scaling relation that also has 
lower scatter than the t raditional M — T re lation is 
the Yx — M500 relation (|Kravtsov et al.l 120061 ). Here, 
M500 is the cluster mass within the radius r§Qo enclos- 
ing an overdensity of 500 relative to the critical density, 
Yx = MgTxi Mg is the cluster gas mass within rsoo, and 
Tx is the mean spectral X-ray temperature of the cluster. 
However, this scaling relation is limited by the assump- 
tions that the gas is both spherically distributed and in 
hydrostatic equilibrium. Our scaling relation offers the 
advantage that it is based on lensing mass estimates that 
are free of these assumptions. 

5. BCG PROPERTIES 

In addition to the interdependencies of many cluster 
properties, properties of the BCG have also been shown 
to correlate with the host cluster. Here we identify the 
BCG in each of our clusters, measure the luminosity of 
each BCG, and examine the correlation between BCG 
luminosity and host cluster mass for our strong lensing 
sample. 

5.1. BCG Determination 

We select each cluster's BCG as the brightest member 
galaxy. Each BCG corresponds to the lens galaxy or one 
of the lens galaxies used to determine the cluster mass 
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TABLE 4 

BCG LUMINOSITIES. 



Cluster 



BCG a 



^K,BCG 
-,11 J.-2 



•t'K.passive.BCG Reference 



C1G 2244-02 

Abell 370 

3C 220.1 

MS 2137.3-2353 

MS 0451.6-0305 

MS 1137.5+6625 

CI 0939+4713 



Gl 



Gl 



1.03 + 0.09 
1.5 ±0.1 

6.4 ±0.4 
8.98 ± 0.09 
4.3 ±0.3 
15 ±2 
1.9 ±0.2 



ZwCl 0024+1652 #362 1.69 ± 0.07 



0.96 ±0.09 
1.3 ±0.1 
5.3 ±0.3 
0.84 ± 0.08 
3.7 ±0.3 
11 ± 1 
1.7 ±0.2 
1.51 ±0.06 



3 

4, 5 



References. — (1 ) Bautz ct al. (1982): (2) Ellingson ct al. (1998)); (3) 
IDe Filippis et al.l 1120031): f4)IKneib et al.l 12003!) ; (5)IMoran et al.l 112005). 



See [Comerford ct al. (2006) for identification of the galaxies by name. 



distribution in IComerford et all (|2006f) . When multiple 
lens galaxies were used to model a single cluster, we iden- 
tify which of the lens galaxies is the BCG in TablelU and 
we also note references that confirm the BCG selection. 

5.2. BCG Luminosity Determination 

For each cluster we have HST imaging taken in 
some combination of the filters F450W, F555W, 
F67 5W, F702W, and F814 W. Using Source Extrac- 
tor (jBertin fc Arnoutsl 119961 ). we measure MAG.AUTO 
magnitudes for the BCG galaxies. We estimate the mag- 
nitude uncertainties by adding in quadrature the error 
in the measured flux and the estimated background sub- 
traction error, which is the product of the area of the 
extraction aperture and the RMS variation of the sub- 
tracted background flux. We calculate the BCG lumi- 
nosities using the available photometry in an observed 
band as the normalization factor on two types of spec- 
tral energy distribution templates, and then compute the 
rest-frame magnitudes and luminosities in several bands 
including if-band. The temp l ates w e use are calculated 
from the Bru zual fe Chariot! (|2003l ) stellar population 
synthesis models with a Salpeter initial mass function. 
The first we use is a fixed-age 10 Gyr old simple stellar 
population, and the second is for a simple stellar popu- 
lation with an age given by an assumed formation red- 
shift of z — 3.0. The latter enables an estimate of the 
passively-evolved BCG luminosity. 

6. THE BCG LUMINOSITY - CLUSTER MASS 
RELATION 

Although it is still unclear how BCGs form, conven- 
tional formation scenarios include galactic cannibalism, 
cooling flows, and mergers during cluster formation (§[T]). 
The evolution of the luminosity of the BCG with the 
mass of the cluster may distinguish between these mod- 
els and offer insight into the formation of BCGs. Semi- 
analytic and numerical simulations of structure forma- 
tion suggest a tight correla tion between BCG luminosi- 
ties and cluster m asses (e.g- JSomervilfe fc Primackll 19991 : 
ICole et al"1l2000l ). and we can parameterize such a cor- 
relation between K-hand BCG luminosities and cluster 
masses M200 by the power law 
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Fig. 4. — The correlation between if-band BCG luminosity and 
cluster mass for our sample of strong lensing clusters. Uncorrected 
luminosities (black points) are fit by the solid line, while lumi- 
nosities corrected for passive evolu tion (open circles) a re fit by the 
dashed line. For comparison, the Lin & Mohr (2004) L — M re- 
lation for the general cluster population is shown as the dotted 
li ne. Our best - fit po wer laws are significantly steeper than that 
of Lin & Mohr (2004), hinting that BCGs in lensing clusters may 
have different formation histories than BCGs in typical clusters. 



Here, we examine the relation between BCG luminos- 
ity and cluster mass for clues about the formation histo- 
ries of BCGs in strong lensing clusters and how their 
formations may differ from the general BCG popula- 
tio n. We represent th e general BCG population with 
the lLin fc Mohr] (|2004D study of 93 BCGs at z < 0.09 in 
the Two Micron All Sky Survey (2MASS). 

For an accurate comparison to the L — M relation 
iLin fc Mohrl (f200"l find fro m 2MASS, we follow their def- 
inition of BCG luminosity. ILin fc Mohrl (|2004h measure 
BCG luminosities in the JT-band using 20 mag arcsec -2 
isophotal elliptical aperture magnitudes for 2MASS, 
called K20 magnitudes. Similarly, we convert to the K- 
band (see § I5.2[) and measure BCG magnitudes using 
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SExt ractor's MAG.AUTO function (jBertin fc Arnoutsl 
119961 ). which has good agreement with 2MASS K20 total 
magnitud es for sources such as BCGs that are bright and 
extended (|Elston et al.ll2006[) . We then convert the mag- 
nitudes into i^-band luminosities as described in § 15.21 
The resultant K-b&nd BCG luminosities, along with the 
luminosities corrected for passive evolution, are given in 
Table |H 

Figure |4] illustrates the correlation between BCG lu- 
minosities and cluster masses M2oo- We find the best-fit 
power law to the data is given by B = 0.97 ± 0.17 and 
f3 = 0.48 ± 0.09 for all strong lensing clusters (solid line 
in FigureH]) and B = 0.93 ± 0.18 and /3 = 0.39 ± 0.10 for 
all strong lensing clusters when the BCG luminosities are 
corrected for passive evolution (dashed line in Figure [4]). 
The similarity of these two results implies that the pas- 
sive evolution of BCG luminosities with redshift has little 
effect on the L — M relation, and more generally there 
is no evidence for evolution in the L — M relation from 
z ~ 1 to z ~ dBrough et al.ll2008h 



For comparison, iLin fc Mohrl pOOl find a best-fit 
power law of B = 4.9 ± 0.2 and /3 = 0.26 ± 0.04 (dot- 
ted line in Figured]), which is consistent with the slopes 
found by analytic estimates and cosmological simulations 
of the growth of central galaxies. Using the galaxy-dark 
matter correlation function to determine host dark mat- 
ter halo masses for ob s ervati onal catalogs of galaxies, 
ICoorav fc Milosavlievid (|2005h find L o c M<° n 3 for halo 



masse s >4x 10 13 h^ 1 M Q . Similarly, IVale fc Ostrikerl 
(2006) determine a correlation of L cx -M^q 8 when they 
combine the subhalo mass distribution derived from sim- 
ulations with an empirical galaxy luminosity function. 
They also find little dependence of the L — M relation 
on waveband. 

From their slope of /3 = 0.26. ILin & Mohrl (|2004f l con- 
clude that while other cluster members may merge with 
BCGs and increase BCG luminosities, such effects are 
not sufficient to fully accou nt for the growth in Lk,bcg 
with cluster mass. Instead, ILin fc Mohrl (|2004f) suggest 
that BCGs must grow mainly through mergers with other 
BCGs brought in when the host galaxy cluster merges 
with other groups or clusters. In addition to the many 
hierarchical structure formation simulations and models 
that support this scenario (g^. iMerrittl [l985t iDubmskl 
119981 : iBovlan-Kolchin et al.ll200"6l ) , there are also observa- 
tions of a pair of ~ L* elliptical galaxies m erging to build 
up th e BCG in a rich cluster at z — 1.26 (Yamada et al. 

Hoo2). 

Our slope /3 is 50% (when luminosities are corrected 
for passive evolution) to 85% (when luminosities are 
not corrected for p assive evolution) steeper than that of 
ILin &: Mohrl i{2()0 1 i. hinting that strong lensing clusters 
may undergo more mergers with groups and clusters, or 
merge with more massive groups and clusters, than the 
average cluster. Both more mergers and mergers with 
more massive systems could account for the initial ev- 
idence for an increase in Lk,bcg with cluster mass we 
find in strong lensing clusters, and would also be con- 
sistent with simulations that suggest strong lensing clus- 
ters are dyn amically more active than the ge neral cluster 
population (|Bartelmann fc S tcinmctzl Il996f) . However, 
the scatter in our L — M relation is significant, and a 
larger sample of strong lensing clusters is necessary to 



draw definitive conclusions about the formation histories 
of strong lensing clusters. 

7. CONCLUSIONS 

We have determined the scaling of cluster mass with 
cluster temperature and the scaling of BCG luminosity 
with cluster mass for eight observed strong lensing galaxy 
clusters imaged with .ffSTand at redshifts 0.3 < z < 0.8. 
We explored cluster concentrations as a means of reduc- 
ing the scatter in the M — T relation and enabling more 
precise constraints on cosmological parameters, and we 
used the L — M relation as an indicator of the formation 
histories of strong lensing BCGs and clusters. Our main 
results are: 

1. The best- fit cluster mass-temperature relation for 
our four dynamically relaxed strong lensing clusters 
is 



E{z) 



M 



2.->()U 



10 14 /ifp M, 



i ■ 



= 1.60 ±3.42 



kT 
5keV 



1.43±1.28 



which is consistent with the theoretical expectation 
of the M — T relation for relaxed clusters, as well as 
the M — T relations determined by other observa- 
tions and simulations. We find that relaxed strong 
lensing clusters do not deviate from the M — T rela- 
tion for the general population of relaxed clusters. 

Significantly, we find an inverse correlation between 
cluster temperature and cluster concentration that, 
when incorporated into the M — T relation, reduces 
the M — T scatter by a factor of 6, from 360% to 
60%. By correcting cluster temperatures accord- 
ing to the temperature-concentration relation, we 
find that the M — T relation given in Equation [7] 
describes not only the relaxed strong lensing clus- 
ters, but the entire cluster population regardless 
of dynamical state. Specifically, the scatter in un- 
relaxed clusters decreases by a factor of 30, from 
500% in the uncorrected M - T relation to 15% 
in the temperature-corrected M — T relation. In- 
corporating concentration effects into the M — T 
relation tightens the M — T relation for all clus- 
ters, which in turn reduces the error in the deter- 
mination of a 8 from cluster counts. Whereas ac- 
curate cluster determinations of a% were previously 
made only with relaxed clusters, concentrations en- 
able the inclusion of unrelaxed clusters. The larger 
cluster samples possible with the inclusion of un- 
relaxed clusters offer yet more precise as estimates 
from cluster observations. 

2. The best-fit relation between BCG luminosity and 
cluster mass for our sample of strong lensing clus- 
ters is 



J K,BCG 



10 11 /if 2 L f . 



0.97 ±0.17 



Mo, 



10 1 



0.48±0.09 



which is ~ 85% steeper than the correlations pre- 
dicted for non-strong-lensing clusters by other ob- 
servations, theory, and cosmological simulations. 
This result supports the current evidence that 
BCGs are built up through mergers with massive 
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galaxies in other groups and clusters, and also hints 
that strong lensing clusters may have more active 
merging histories than typical clusters. A larger 
sample of strong lensing clusters might enable more 
definite conclusions about the formation histories 
of strong lensing clusters. 

Accurate cluster mass measurements and full use of 
the range of cluster property interdependencies are key 
components in the calibration of clusters as tracers of cos- 
mological parameters. As we have shown, gravitational 
lensing enables the most direct measurements of cluster 
mass, without assumptions about the cluster's dynamical 
state that are inherent in other methods. We have also 
shown that the correlation between cluster temperature 
and concentration can significantly reduce the scatter in 



the cluster M — T relation, enabling more precise esti- 
mates of a%. It may be that other cluster scalings can 
be effectively combined to reduce the error on additional 
cosmological parameter estimates. 
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